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ABSTRACT 

We present the results on period search and modehng of the cool DAV star KUV 
02464+3239. Our observations resolved the multiperiodic pulsational behaviour of 
the star. In agreement with its position near the red edge of the DAV instability 
strip, it shows large amplitude, long period pulsation modes, and has a strongly non- 
sinusoidal light curve. We determined 6 frequencies as normal modes and revealed 
remarkable short-term amplitude variations. A rigorous test was performed for the 
possible source of amplitude variation: beating of modes, effect of noise, unresolved 
frequencies or rotational triplets. Among the best-fit models resulting from a grid 
search, we selected 3 that gave I = 1 solutions for the largest amplitude modes. These 
models had masses of 0.645, 0.650 and 0.680 Mq. The 3 'favoured' models have Mh 
between 2.5xl0~'^ — 6.3xlO~®Af* and give 14.2 - 14.8 mas seismological parallax. 
The 0.645 Mq (11 400 K) model also matches the spectroscopic \ogg and Toff within 
1 (7. We investigated the possibility of mode trapping and concluded that while it can 
explain high amplitude modes, it is not required. 
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1 INTRODUCTION 

Since Arlo Landolt's (1968) discovery of the first DAV 
(HL Tau 76) we have learned much about pulsating white 
dwarfs. They are otherwise normal white dwarf stars, and 
they represent the fi nal stage of about 98% o f the stars (see 
the recent review of IWinget fc Kepleill200^ and references 
therein), so it is a key issue to gain insight into their interi- 
ors. While white dwarfs are relatively simple to model, some 
observed phenomena still challenge theory. 

One unsolved mystery is the variability of the pulsa- 
tional modes' amplitudes and frequenc ies of PNNV, cool 
DBV and cool DAV stars (|Handleil2003al '). While we can find 
very stable modes in certain white dwarfs which can be used 
to measure evolutionary effects (cooling and/or contraction, 
e.g. hot DAVs G117-B15A and R548 - 'The Most Stable Op- 
tical Clocks Known': Mukadam, Kepler & Winget 2001 and 
references therein), in some cases short-term variabilities 
make it difficult to determine even the pulsational modes. 
An example for this latter is the case of PC 1115-1-158 (DBV) 



E-mail: bognar[at]konkoly[dot]hu 



with i ts remarkably unstable amplitude spectra (|Handleij 
l2003bl ). 

In some cases we encounter strange, sudden, short-term 
variations in the pulsational behaviours of stars: for exam- 
ple the 'sforzando' effect observed in the DBV star CD 358 
in 1996, where the nonlinear light curve of the star changed 
into a high ampl itude, remarkably s i nusoidal one for a shor t 
period of time (|Kepler et al.1 120031 . iProvencal et al.l bOOgl ). 
In the case of PG 1456-1-103 its luminosity variations al- 
most d isappeared just be fore the Whole Earth Telescope 
(WET , [Na"ther et al.|[l990l ) run started on the star (|HandIei) 
l2003al ). The frequency and amplitude variations of multiplet 
components (mostly fo r high k modes) observed in GD 358 
IProvencal et al.ll2009l ) are further examples for short-term 
variations. 

Possible explanations for these observed phenomena 
could be non- linear mode coupling, convection/pulsation in- 
teractions, mode excitation, and beating of different pulsa- 
tion modes - such as unresolved rotational splitting. To find 
the right explanation(s) for a certain case is a great chal- 
lenge. A crucial first step is to identify the periods of the 
true pulsation modes of a white dwarf. 

One possible way of mode identification when strong 
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amplitude variations occur is to observe the star in many 
seasons, determine all the observed periods from the differ- 
ent runs and search for the (nearly) equidistant period spac- 
ings among them, which can be con secutive overtones with 
the same horizontal degree (Z) value. iKleinman et alj (|f 9981 ) 
used this method for the first time for a DAV star (G29-38) 
where recurring sets of normal modes were identified. 

Since KUV 02464+3239 was found to be a luminosity 
varia ble DA star near the red edge of the DAV instability 
strip (|Fontaine et al.|[200lf ) , we expect to find the short-term 
variabilities which characterize the pulsation of similar stars. 
Using a 2900 s-long light curve. [Fontaine et al.l (|200fl ') deter- 
mined a quasi-period of ~832 s for the star's light variation. 
The long period, large amplitude variations they found with 
strongly non-sinusoidal light curve along with the presence 
of harmonic peaks in the frequency spectrum are consistent 
with the behaviour of a cool DAV star. However, up to now 
only the results of this short discovery run had been pub- 
lished about this target. 

In this paper we present the results of observational 
analysis of KUV 02464-1-3239. We describe our observations 
and data reduction process in Sect. [2] The results of Fourier 
analyses and frequency determination tests can be found in 
Sect. |3] Sect. |4] contains our tests on amplitude variations. 
We discuss the results of asteroseismological investigations 
(modeling, results for stellar parameters) in Sect.[5]and give 
a brief summary in Sect. (6] 



2 OBSERVATIONS AND DATA REDUCTION 

We observed KUV 02464-1-3239 on 20 nights between Oc- 
tober 2006 and February 2007. The observations were 
made with a Princeton Instruments VersArray;I300B back- 
illuminated CCD camercQ attached to the Im RCC tele- 
scope at Piszkesteto mountain station of Konkoly Observa- 
tory. The log of observations is in Table [T] We did not use 
any filter ('white light observations') to maximize the pho- 
ton count from this relatively faint target (rrii, = 16™07). 
The back-illuminated detector we used is more sensitive to 
bluer wavelengths than a front-illuminated one. Note that 
the amplitude values of the Fourier light curve solution cor- 
responding to these white light measurements differ from 
filtered ones. 

Standard IRAE0 routines were used during the reduc- 
tion process; the frames were bias, dark and fiat corrected 
before we performed aperture photometry with the IRAF 
DAOPHOT package. The observational times of every data 
points were calculated to Barycentric Julian Date (BJD). 

Different apertures were used for every night. For a 
given night, we chose the aperture size that resulted in the 
lowest scatter for differential light curves of constant stars 
in the field. 

We used 10 or 30 s exposure times depending on the 



^ http:/ /www. princctoninstrumcnts.com/products/imcam/ 
versarray / dsheet .aspx 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Table 1. Journal of observations of KUV 02464-1-3239. Subset 
No. refers to subsets in Sect. 13.21 A' denotes the number of points 
remained after averaging (see details in text) and ST means the 
time span of the runs including gaps. 



Run Subset Date Start time N &T 

No. No. (UT) (BJD-2 450 000) (h) 



01 
02 
03 
04 
05 
06 
07 
08 
09 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



Total: 8227 97.80 
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Figure 1. Variable and comparison stars in the CCD field. The 
three brightest stars (CI, C2, C3) were selected as a reference 
system for the differential photometry. 

weather conditions during the observing runs. For 10 s ex- 
posure we averaged the data into 30 s bins to keep all the 
data with the same bin, although, we are aware of that for 
multiperiodic pulsators with nearby pulsations it does intro- 
duce an extra uncertainty. 

We performed aperture photometry on the same 22 
stars for every night. Possible comparison stars were checked 
for variability and longer-term trends (caused by colour ef- 
fects of stars with different spectral types) in an iterative 
way. This process resulted in our choosing the three bright- 
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Figure 3. Normalized difTcrcntial light curves of the 20 runs obtained on KUV 02464+3239 between October 2006 and February 2007. 
Time is given as B JD-2 450 000. 




Figure 2. Differential light curve of the check star (C3) to the 
average of the reference stars obtained on JD 2 454 086. 



est stars marked in Fig. [T] as a reference system for the dif- 
ferential photometry. The quality of this method is shown 
by the differential hght curve of a check star (C3) compared 
to the average of the reference stars in Fig. (2] 

Following the traditional method of photoelectric pho- 
tometry, we obtained first-order extinction coefficients for 
the nights which were long enough. The values lay in the 
range of 0.07-0.39 in white light. 

The colour term of extinction was checked. However, we 
could not precisely determine the second order extinction 
coefficient for each night. In addition to the iterative pro- 
cess, we obtained BVRI photometry to estimate the spec- 
tral types of the comparison stars selected in the iterative 
process. The relative colour indices were checked. The dif- 
ferences between the B-V colour indices are higher than 



~0.5 which shows immediately that the comparison stars 
have much later spectral types than KUV 02464-1-3239. Nat- 
urally, this deficiency often resulted in trends in the differen- 
tial hght curves and accordingly significant signals in the low 
frequency region of their Fourier Transform (FT) . Although 
the low frequency signals do not exceed 1-6 mmag, to get a 
homogeneous dataset for the variable star, polynomial fits 
were performed in the last step of the reduction process. It 
is widely applied in white dwarf research using multichannel 
photometry. The light curves obtained for each run can be 
seen in Fig. O 



3 FOURIER ANALYSIS 

Since the star was not well studied before, our first aim was 
to find the frequency content. Analyses were made by use of 
the MuFrAn (Multi-Frequency Analyzer) package (|Kollathl 
Il99d : ICsubrv fc KoUathl 12001 ) . MuFrAn provides efiicient 
tools for frequency determination with its standard ana- 
lyzer applications (FFT, DFT, linear and nonlinear fitting 
options) and graphic display routines. The program handles 
unequally spaced and gapped observational data. Frequency 
values are given in cycle/day (c/d) units. 
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Table 2. Frequencies detected with S/N ^4 values in the Fourier 
Transforms of subsets. Period and amplitude values are given for 
completeness. 
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Figure 4. Amplitude spectra of the subsets. Numbers of subsets 
are indicated in the left-upper corners of each panels. Window 
functions are given in the inserts. Remarkable amplitude varia- 
tions can be seen from one interval to another. 



We followed the standard steps of a pre-whitening pro- 
cess to get frequency, amplitude and phase values utilising 
the light curve's Fourier spectrum. In pre-whitening pro- 
cesses the question is always when to stop iterating. In decid- 
ing whether a peak belongs to a real pulsation frequency, we 
used the widely accepted approximation: if a peak reaches 
the threshold of signal-to-noise r atio of ~4 it has a small 
probability of being due to noise (|Breger et al.lll993l ). 

The tools of the time string analysis program Period04 
l|Lenz fc Bregeill2005l ) were used to calculate S/N of the in- 
dividual peaks. The program determines the signal (S) value 
as the amplitude of the peak after the least-squares fit. The 
noise (N) is given as an average amplitude calculated in 
a frequency range that encloses the peak in the residual 
spectrum. We determined S/N values from the ±25 c/d fre- 
quency range of peaks after pre-whitening the original FT. 



3.1 Investigation of individual nights 

Because the pulsation periods are short compared to the 
typical run length, the one-night long observing runs provide 
a long enough time base to investigate the basic pulsational 
characteristic of the star, even in a case if not all the modes 
are resolved. In addition, a night by night analysis allows us 
to follow the changes in pulsation amplitudes. 

Considering the threshold of S/N~4 for the peaks, 
our results show that the pulsation frequencies of KUV 
02464+3239 are at ~87, 97, 104, 111 and 139 c/d. In one case 
we also find a very significant (S/N=5.4) signal at ~70c/d. 
During the pre-whitening process we did not take into ac- 
count peaks under 60 c/d in the FTs. This limit is an overes- 
timation of the frequency range filtered by polynomial fitting 



Frequency 


Period 


Ampl. 


b/iN 


(c/d) 




(mmag) 




Subset 1 








90.660±0.004 


953.01 


12.8 


4.9 


99.251±0.003 


870.52 


17.4 


6.6 


103.790±0.002 


832.45 


30.6 


11.9 


114.103±0.004 


757.21 


10.0 


4.0 


Subset 2 








99.567±0.028 


867.76 


23.6 


7.1 


104.880±0.036 


823.80 


20.2 


6.5 


111.452±0.045 


775.22 


15.1 


5.6 


139.227±0.063 


620.57 


11.6 


6.0 


Subset 3 








69.475±0.047 


1243.61 


14.7 


5.4 


87.264±0.044 


990.09 


16.3 


6.0 


98.025±0.053 


881.41 


16.1 


6.0 


104.358±0.047 


827.92 


18.0 


7.3 


Subset 4 








*86.579±0.004 


997.93 


9.5 


5.5 


86.996±0.002 


993.15 


18.4 


10.6 


*89.142±0.003 


969.24 
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6.2 


*90.100±0.001 


958.94 


15.2 


8.8 


98.797±0.002 


874.52 


9.1 


5.5 


104.036±0.002 


830.48 


10.2 


6.5 


139.954±0.005 


617.35 


3.8 


4.7 


173.960±0.004 


496.66 


4.7 


4.3 



* Closely spaced frequencies to the dominant peak. Signals at 
89.14 and 90.10c/d are close to the Ic/d aliases of each others. 
The peak at 90.1 c/d has an unreliably high amplitude solution 
only when we fit at least the first six frequencies of this subset. 

and this means that any possible pulsation frequency below 
this level remained undiscovered. 

Our analysis revealed that remarkable amplitude varia- 
tions can happen from one night to anothe r as in many cas es 
of PNNV, cool DBV and cool DAV stars (|Handleij|20o"3al ). 



3.2 Analysis of reference intervals 

The dataset can be grouped into four subsets of selected 
nights. The longer time base means better spectral resolu- 
tion and with the proper selection of nights, the time base 
still remains short enough not to obscure the possible short- 
term variabilities. 

The first four nights (JD 2 454 014-2 454 019) were se- 
lected from October (subset No. 1). We got the longest time 
string during the eighth run (JD 2 454 068) so we used this 
run as a reference for November (subset No. 2). In December 
we had the chance to observe KUV 02464+3239 in two con- 
secutive weeks. In view of the results of individual nights' 
Fourier analyses we selected two intervals; the first is the 
run on JD 2 454 081 (subset No. 3), the second consists of 
four nights (subset No. 4; JD 2 454 084-2 454 089, except 
the short run on JD 2 454 087). The latter runs have similar 
FTs but the run on JD 2 454 081 has a different FT. 

Fig. |4] shows the FTs of our subsets. The variations in 
amplitudes of peaks in the FTs are obvious. Frequency val- 
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Table 3. Frequency and amplitude values of the 6 accepted pul- 
sation frequencies for each subset. Amplitudes were determined 
by fitting with these frequency components only. Subset No. refers 
to subsets of Sect. 13.21 



Subset 






Frequency 




No. 








(c/d) 




1 






99.25 


103.79 




2 






99.57 


104.88 111.45 


139.23 


3 


69.48 


87.26 


98.02 


104.36 




4 




86.99 


98.80 


104.04 




Amplitude (mmag) 


1 






17.28 


29.56 




2 






23.61 


20.24 15.12 


11.59 


3 


14.74 


16.25 


16.07 


17.96 




4 




20.83 


8.83 


10.14 





ues determined after pre- whitening are given in Table[2] S/N 
values of peaks were also calculated. Noise levels change from 
lower to higher frequencies between 2.6-2.5 mmag (subset 1), 
3.3-1.9 mmag (subset 2), 2.7-2.5 mmag (subset 3) and 1.7- 
Immag (subset 4). Standard deviations of frequency values 
are determined by Monte Carlo (MC) simulations, getting 
solutions for each frequency in each subset. We generated 
synthetic hght curves (by adding Gaussian random noise) 
corresponding to the frequencies can be seen in Table [21 
After the non-linear least squares fitting of 100 synthetic 
datasets of each subset we determined the standard devia- 
tions of the frequencies. 



3.2.1 Frequency determination test 

To support our final frequency determination we performed 
a robustness test. Synthetic hght curves were generated us- 
ing frequencies, amplitudes and phases of the hght curve 
solutions with Gaussian random noise added. The applied 
noise levels were scaled from 1 to 4 times the residual scat- 
ters (a) of the light curve solutions. Independent frequency 
analyses were carried out for each dataset (1) to check how 
accurately we can get back the input frequencies from the 
noisy datasets and (2) whether we can identify the same 
frequencies that were found to be significant in the original 
datasets. 

Testing our first subset we faced a ±1-2 c/d alias prob- 
lem as we raised the noise level, especially in the case of 
frequencies with lower amplitudes. Frequencies were found 
rather far from the input values. At the 4 a level we could 
identify 2 frequencies (99.25, 103.79 c/d) out of the 4 found 
to be significant earlier. In the case of our second and third 
subsets, the amplitudes were rather high and we found all 
4 of the frequencies identified before. Adding noise to the 
fourth subset caused alias problems from the 2.5 a level up- 
ward and the lowest amplitude peaks disappeared. We could 
find only 4 frequencies (86.99, 90.10, 98.79 and 104.03 c/d) 
within the most noisy time string. 

Table[3]shows the results of our robustness test with the 
finally accepted frequencies. Peaks at ^104 and 99 c/d were 
found in all subsets. Other four frequencies (at ~69, 87, 111 
and 139 c/d) were rather high amplitude peaks in one or two 
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Figure 5. Frequecies out of the 6 accepted ones found at least 
in two subsets. We calculated the width of error bars by the fre- 
quency range within the peaks were found during our robustness 
test (see Sect. 1372. II for details). 



segments. These 6 frequencies were determined unambigu- 
ously by the analyses of subsets and perhaps characterize 
the pulsation of KUV 02464-1-3239 within some amplitude 
limit. 

As can be seen in Fig. [4] and in Table [3l remarkable 
amplitude variations occur from one subset to another. A 
signal at ^87 c/d became dominant in December while the 
amplitude of the peak at 104 c/d decreased by 66%. The 
peak at 99 c/d first increased its amplitude by 37% then 
decreased by 63%. 

We represent the frequencies which were found at least 
in two subsets together in Fig. [5] The groups of frequen- 
cies are well-separated, there is no problem with overlap- 
ping sidelobes. The error bars give the uncertainty of our 
findings on a certain value in a subset. Since they overlap, 
it is clear that the two peaks around 87 c/d can have differ- 
ent values because of the uncertainties. The error bars also 
overlap in the case of peaks at ~104c/d. The four peaks 
around 99 c/d represent an interesting case. Since the three 
well-determined peaks are more widely separated than the 
error bars, we can say that they are distinct. The reasons 
could be a frequency change on a short time scale (from one 
subset to the other - within one month). Another expla- 
nation is that they are independent modes that are always 
excited and their amplitudes change from subset to subset. 
Alternatively, they might be components of an unresolved 
rotationally split mode. A comparative analysis on a longer 
time base could clear up the real situation (if it is a regu- 
lar behaviour). At this early stage of the interpretation of 
KUV 02464+3239 we mention this behaviour but we give a 
general solution based on the whole dataset. 

3.3 Analysis of the whole dataset 

The long time base of the whole dataset gives more precise 
frequency values if we assume the frequencies are stable over 
the time base. In this case the short-term variations in ampli- 
tudes are obscured but the frequencies are well-determined. 

The successive pre-whitening steps of the whole dataset 
can be seen in the panels of Fig.[Sl Table 3] shows the values 
of frequencies determined from the whole dataset. Standard 
deviations are calculated by MC simulations. Peaks at /i - 
/e are marked in the first and second panels of Fig. [S] (upper 
segment of Table [4]) . They represent the 6 frequencies that 
come out of the final analysis. The residual spectrum in the 
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Figure 6. Frequency analysis of the wliole dataset. Tlie panels 
show the successive pre- whitening steps, the window function is 
given in the insert. The 6 finally accepted frequencies are marked 
in the first and second panels. The third panel shows the residual 
spectrum with definite signals remained. The spectrum after prc- 
whitening with 13 frequencies can be seen on the last panel. It is 
still not white noise. 

third panel is definitely not white noise. With further anal- 
ysis, we found 7 other frequencies (fr-fia) in the residual 
spectrum and present them in the lower segment of Tabled 
The last panel of Fig. |6] shows the FT after pre- whitening 
with 13 frequencies. 

Amongst fy, fg and fg, additional pulsation modes could 
be present, but we cannot say which ones are real modes, fr 
is near to the -1-1 c/d alias of the /2 peak obtained in subset 
3. There is no sign of fg in any subsets, fg can be a remainder 
of the group of four peaks near /s obtained by the subsets. 
The two higher frequencies in Table [3] (99.25, 99.57 c/d) ap- 
pear as /a and the two lower frequencies (98.02, 98.80 c/d) 
can produce fg in the whole dataset. The double structure 
around 139 c/d might be a sign of rotational splitting: 5 f 
= 0.5 c/d = 5.79 /iHz. Assuming that these peaks belong 
to high-overtone (k ^ 1) I = 1, m = —1,0 or m = 0, 1 
modes the rotation period of the star would be ~ld. Con- 
sidering the peaks at 104 c/d, their frequency separation is 
only 5f = 0.13 c/d = 1.5 ^Hz which corresponds to a ~3.8d 
rotation period under the assumptions above. However, re- 
garding the large amplitude of the closely spaced peak, fio 
could also be an independently excited mode with another 
I value. Determination of the frequency at 111.09 c/d (fn) 
was ambiguous because of aliasing. The peak at 173.96 c/d 
is consistent with being the first harmonic of the dominant 
mode. 

For asteroseismological investigations, we use only the 
6-frequency solution supported by the analyses of subsets. 
We consider these frequencies as a set of normal modes. 
A dataset with better coverage and higher signal-to-noise 



Table 4. Frequency, period and amplitude values of 13 frequen- 
cies derived by the whole dataset. Frequencies at /i - /g arc ac- 
cepted as normal modes by the analyses and test of subsets. Fre- 
quency and amplitude values of /i - /g were determined by fit- 
ting with these components only. Frequencies and amplitudes of 
fy^fia were calculated by fitting with all the 13 frequencies. Am- 
plitude values of /i — /e for the 13 frequency solution are given in 
parentheses. 



Frequency Period Ampl. 

(c/d) (mHz) (s) (mmag) 



/i 69.1060 ±0.0003 799.838 1250.253 4.4 (4.4) 

/2 86.9879 ±0.0001 1006.804 993.242 13.2 (14.4) 

h 99.7516 ±0.0001 1154.532 866.151 9.5 (11.5) 

/4 104.2617 ±0.0001 1206.733 828.684 11.6 (11.4) 

/s 111.1095 ±0.0002 1285.989 777.611 5.5 (8.0) 

/6 139.5157 ±0.0003 1614.765 619.285 4.0 (3.9) 



h 88.1611 1020.383 980.024 8.4 

fg 89.7100 1038.311 963.103 7.3 

fg 98.8533 1144.136 874.022 6.5 

fio 104.3915 1208.235 827.654 10.2 

fii 111.0865 1285.724 777.772 5.8 

fi2 139.0162 1608.983 621.510 3.7 

fi3 173.9642 2013.474 496.654 3.2 



ratio would enable us to determine additional real pulsation 
modes. 



4 TESTS ON AMPLITUDE VARIATIONS 

The Fourier analysis method presumes stationary ampli- 
tude and frequency values while the behaviour of KUV 
02464+3239 seems to be different. We do not know why the 
observed amplitudes of the star vary so much. To provide 
constraints for possible explanations, we performed tests for 
the amplitude variability. 

We generated synthetic hght curves corresponding to 
the whole dataset using the frequency, amplitude and phase 
values of either the 6 accepted frequencies or with the addi- 
tional closely spaced frequencies around the 87 c/d dominant 
mode. In every case, we analysed the four subsets by nonlin- 
ear least-squares fitting using the original frequency values 
as initial parameters. 

As our first test (a) we generated a synthetic hght curve 
using only 6 frequencies and we did not add noise to the 
dataset. The analyses of subsets did not show noticeable 
amplitude variations from one subset to another. This result 
shows that beating of the 6 modes cannot explain the ap- 
parent amplitude variations. In our second test (fe) we added 
Gaussian random noise to the synthetic hght curve using the 
residual scatter of the whole light curve's 13-frequency solu- 
tion. Our analysis showed that noise had a slight infiuence on 
the amplitudes (mainly for the lower amplitude modes) and 
cannot be responsible for the large observed variability. In 
our subsequent tests (c - h) we generated the synthetic hght 
curves using 6+2 frequencies where the two additional fre- 
quencies were close to the dominant mode {5f = ±0.185 c/d 
= 2.14 /^Hz). This separation was selected to be close to the 
resolution limit of the longest subset (subset 4) and this 
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Table 5. Results of tests on amplitude variability of the 87c/d 
mode. At test cases a - h we give relative amplitude differences 
with respect to the synthetic light curves' input value. Sec details 
in Sect, g] 



Subset 



Test cases 
c d e f g 

Amplitude variation (%) 



1 





+4 


+28 


+9 


-12 


-3 


+10 


-3 


2 





-6 


+60 


+37 


-52 


+11 


+1 


-9 


3 





+15 


-87 


-64 


-17 


-7 


-28 


+6 


4 





-5 


+22 


+11 


-35 


+7 


+14 


-4 



could be a reasonable value for rotational splitting. We set 
the amplitudes of the triplet's side components to be the half 
(c - e) and one tenth (/ - h) of the 87c/d mode's amplitude. 
The phase differences were ±22.5° (c, /), ±45° (d, g) and 
±90° (e, h) with respect to the central peak's phase at the 
initial epoch (JD 2 454 014). With these tests we simulated 
the effects on amplitudes of closely spaced frequencies which 
might remain unresolved in the subsets. Such peaks could 
be detected by the analysis of the whole hght curve if they 
had large enough amplitudes so as not to vanish in the other 
peaks and the noise. 

Table [5] summarizes the test results on the ampli- 
tude variations of the dominant peak. In some cases we 
got large amplitude variations: disappearance and almost- 
disappearance of the peak in subset 3 (tests c and d). In 
test c we fitted a peak at 86.3 c/d (the input value was 
86.99 c/d) with 1.5mmag amplitude. This means an ampli- 
tude decrease by 87% if we accept the 86.3 c/d peak as the 
representative of the original one. In test d we got a 64% de- 
crease in amplitude. Triplet side components with smaller 
amplitudes caused smaller, but noticeable variations. The 
maximum was a 28% decrease in the amplitude of the dom- 
inant peak (test g). 

Amplitude variations were detected for other modes as 
well. The presence of noise and beating with triplet compo- 
nents can be the explanation since we generated the syn- 
thetic hght curves using constant amplitude values. The av- 
erage of the variations of the 5 modes is 15%. 

We conclude that the observed large amplitude variabil- 
ity can be simulated with a triplet assuming certain phase 
relations and relatively high amplitudes for the side com- 
ponents. However, these special phase relations and high 
amplitudes of side lobes to the main peak are rather im- 
probable. A real amplitude variation seems to be plausible. 



5 ASTEROSEISMOLOGY OF KUV 02464+3239 

Our main goal in asteroseismology is to find stellar mod- 
els that match the observed properties of the star with the 
highest precision. The most important observed values are 
the periods of hght curves. At first we are looking for mod- 
els which produce pulsation periods with the lowest dif- 
ferences to the observed ones. We used the White Dwarf 
Evolution Code (WDEC) orig inally written by Martin 
Schwarzschild and modified by iKutter fc Savedofj (|l969l ). 



Lamb fc van HornI (Il975l ) IWingetl (Il98ll). iKawaleil (Il986h . 
WoodI (I1990I). iBradlevI (Il993l). iMoiitgomeryl (|l998l ) and 

Bischoff-Kim. Montgomery fc Winged ( 2008al ) to generate 



equilibrium white dwarf models each with a given set of 
input parameters. For each model, we computed a set of pe- 
riods based on the adiabatic e quations of non-radial stellar 
oscillations (|Unno et al.lll989l ). We computed grids of such 
models, varying the input parameters. 

In WDEC we start with a hot (~100 000K) polytrope 
model. The code evolves it down to the requested temper- 
ature and the model we finally get are thermally relaxed 
solution to the stellar structure equations. To get pulsation 
periods of a model we used the integrated form of the WDEC 
given by Metcalfe (2001), which includes a pulsation code. 

W e used the equation-of-state (EOS) tables of 
LambI (Il974l') in the core and the EOS tables of 



Saumon. Chabrier fc van HornI (|l995l ) in t he envelope of the 



star. W e used OPAL opacities updated bv llglesias fc Rogers! 
(| 19961 ). The convection was treated within the Mixini 



Lengt h Theory (MLT) according to iBohm fc Cassinell 
(|l97ll ) and with the a = 0.6 paramctrization suggested by 
model atmosphere calculations of Be rgeron ct a l. (1995tJ). 
Precise value of a has very little effect on the pulsation 
periods of models (Bischoff-Kim, Montgomery fc Winget 
2008b). Treatment of the hydrogen/ helium transition zone 
was done by equilibrium diffusion calculations contrary to 
the treatment of the helium/carbon transition layer which 
was parametrized. 



5.1 The model grid 

As we discussed in Sect. 13.31 we used the 6-frequency 
light curve solution for the seismological analysis of KUV 
02464+3239. Table g] shows that we found modes only 
in the long period regim e (between 619s and 1250s). 
iBischoff-Kim et al.l (|2008bl ) shows that low order modes are 
especially sensitive to the mass of the hydrogen layer. 

The best-fitting models have to fulfil two criteria: (1) 
to give back the period values closest to our observed values 
and (2) to provide physically acceptable solutions for effec- 
tive temperature and logg. Finding solutions according to 
the first criterion we used the f itper tool developed by Kim 
(2007). The program calculates the r.m.s. values using the 
following equation: 



:>obs\2 



iV 



(1) 



where denotes the number of observed periods. 

We built our model grid varying 5 input parameters 
of the WDEC: Tcs, M, , Mh, Xo (central oxygen abun- 
dance) and (the fractional mass point where the oxy- 
gen abundance starts dropping). For our first scan we fixed 
the mass of the helium layer at Mne = 10~^Mj,. We built 
our grid and searched for acceptable solutions in the pa- 
rame ter space determined by spectroscopic results on the 
star. iFontaine et~al] (|200ll ) derived Toff = 11 290 K and 
log 5 = 8.08 values for KUV 02464+3239 with the esti- 
mated fl ux calibration uncert ainties of ±200 K and 0.05 dex 
in log g (|Fontaine et al.ir2003l ). However, this uncertainty es- 
timate does not include contributions from different spec- 
tra, modeling uncertainties, and different fitting procedures. 
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Table 6. Best-fitting sol utions for KUV 02464 +3239 in the mass range 0.525 - 0.74 Mq. We represen t the r esults of former spectroscopic 
observations (Tcff , logo. [Fbntaine et al.|[200lh . the expected seismological mass derived by iBradlevI l|l996l ) and the observed pulsation 
periods as reference values. Models that have masses within 1 cr in logg range are typeset in boldface. 



A/,/A/q, (logg) T|,fj- (K) -logA'/jj Xa -X^fm Model periods in seconds {l,k) Cj-.m.s. (s) 



0.525, (7.85) 


11400 


6.9 


0.7 


0.5 


620.4 

(1,8) 


777.1 

(1,11) 


830.3 
(1,12) 


865.4 
(2,24) 


992.3 

(1,15) 


1250.7 

(1,20) 


0.95 


0.525, (7.85) 


11600 


7.2 


0.5 


0.4 


620.3 
(2,16) 


776.3 
(1,11) 


830.7 
(1,12) 


866.0 
(2,24) 


992.6 
(1,15) 


1251.0 
(1,20) 


1.12 


0.575, (7.97) 


11000 


6.2 


0.7 


0.5 


619.8 
(2,17) 


777.4 
(1,12) 


831.0 
(1,13) 


868.7 
(2,25) 


994.0 
(1,16) 


1250.5 
(1,21) 


1.45 


0.585, (7.99) 


11000 


6.4 


0.5 


0.4 


617.9 

(2,17) 


775.4 

(1,12) 


827.1 

(1,13) 


866.9 

(2,25) 


991.6 
(1,16) 


1250.5 
(1,21) 


1.44 


0.615, (8.03) 


11800 


4.0 


0.7 


0.3 


616.4 
(1,12) 


777.4 
(1,16) 


828.8 
(2,31) 


865.6 
(1,18) 


991.1 
(1,21) 


1249.7 

(2,48) 


1.51 


0.625, (8.04) 


11000 


7.4 


0.5 


0.2 


620.6 
(1,9) 


780.0 

(1,12) 


827.4 

(2,24) 


866.8 
(1,14) 


994.1 
(2,29) 


1252.1 
(1,21) 


1.50 


0.645, (8.07) 


11400 


5.2 


0.9 


0.2 


618.1 
(2,21) 


779.7 
(1,15) 


828.2 
(1,16) 


865.7 
(1,17) 


991.9 
(1,20) 


1251.8 

(2,45) 


1.33 


0.650, (8.08) 


11800 


4.6 


0.6 


0.1 


620.3 
(1,12) 


776.6 
(2,29) 


827.1 
(1,17) 


866.2 
(1,18) 


992.8 
(1,21) 


1249.6 
(2,48) 


0.93 


0.680, (8.13) 


11800 


5.0 


0.5 


0.1 


618.8 
(1,12) 


778.1 

(2,29) 


826.6 
(1,17) 


865.2 
(1,18) 


992.1 
(1,21) 


1251.9 

(2,48) 


1.26 


0.685, (8.14) 


11400 


4.8 


0.6 


0.4 


618.7 
(1,12) 


778.6 
(1,16) 


826.4 
(1,17) 


866.4 
(2,32) 


994.1 
(2,37) 


1249.9 
(2,47) 


1.12 


0.725, (8.2) 


11800 


5.8 


0.7 


0.4 


620.4 
(1,12) 


778.5 

(1,16) 


828.4 
(1,17) 


865.6 
(1,18) 


995.3 
(1,21) 


1249.5 

(1,27) 


1.11 


0.725, (8.2) 


11200 


5.4 


0.8 


0.4 


616.9 
(1,12) 


778.5 
(1,16) 


828.9 
(1,17) 


865.2 
(1,18) 


992.7 
(2,37) 


1251.6 
(2,47) 


1.26 


0.740, (8.23) 


11800 


6.0 


0.6 


0.4 


620.0 
(1,12) 


778.9 

(1,16) 


828.1 
(1,17) 


867.0 
(1,18) 


991.3 

(1,21) 


1250.7 
(1,27) 


1.09 


Reference values: 
























0.65, (8.08) 


11290 








619.3 


777.6 


828.7 


866.2 


993.2 


1250.3 





We make some allowance for this by covering a range of 
±500 K in Tcft and ±0.1 dex in logg. According to this ef- 
fective temperature value, our grid covers the range 10 800 
to 11 800 K (~11290K±2.5o-). Estimating the mass range 
we have to cover, we used the se ismological masses deter- 
mined for DAs by iBradlevI (|l996l ). According to his results 
the mass of a DA white dwarf with logi; = 8.08 is about 
0.65 Mq. Our grid covers the M, = 0.525 - 0.74 A/q mass 
range (loggr =~ 7.9 — 8.2) which means that we investi- 
gated model solutions at \ogg = 8.08±2(t. We varied the 
mass of the hydrogen layer between lO"** — 10~*A^, . The 
core parameters were changed between Xo = 0.5 — 0.9 and 
Xim = 0.1 — 0.5 taking into a c count the carbon-oxygen pro- 
files derived by ISalaris et al.l (|l997l ) . Step sizes were 200 K 
(Teff), 0.005 Af0 (Af,), 10-°'^Af. {Mn) and 0.1 (Xo and 



5.2 Best-fit models 

While running the WDEC we allowed the modes to be Z = 1 
or 2. Since there are no previous asteroseismological analy- 
ses on this star or other contrainst on the I values of modes 
we let all 6 modes be Z = 1 or 2 for the fitting procedure. 



The variations in pulsation amplitudes rule out a mean- 
ingful application of differe nt weights on the amplitudes 
ICastanheira fc Keple3l2008h . 

Assuming better visibility of Z = 1 mo des over I = 2 
modes (see e.g. ICastanheira fc Kepleij [ioOSi and references 
therein), we prefer models that gave at least 3 / = 1 solu- 
tions with low Or.m.s. values. We summarize the parameters 
of the best-fitting models in Table (6] It shows the model 
parameters, calculated pulsation periods with their I and k 
values and the corresponding ar.m.s. value. We find that the 
mode identifications of the best-fitting models are typically 
quite different from model to model. This is due to the fact 
that over the span of 400 K, the periods of the relatively high 
overtone modes observed can change by about 40 s. As a re- 
sult, the cooling of models for a given structure causes differ- 
ent ? = 1 or 2 modes to fit a given observed period best. Note 
that since our current models use abo ut 400 zones, the un- 
certa inty of the model periods is ~1 s (jSrassard fc Fontaine! 
l2005h . therefore, differences in periods below 1 s are within 
the 'noise' of our modeling accuracy. The presence of shorter 
period modes could give stronger constraints on the stellar 
structure, resulting in fewer solutions with low Or.m.B.- How- 
ever, within the '^l s uncertainty of our modeling we can 
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Figure 7. The 13 models of Tablc|6]in the T<,ff - A/, plane. Closed 
and open circles denote the models within and out of the 1 a log g 
range determined by spectroscopy. The black square denotes the 
spectroscopic solution with its uncertainties given by Fontaine et 
al. (2001, 2003). The grid with dashed lines corresponds to the 
step sizes of our model grid. 



map the possible groups of models in the parameter space 
determined by spectroscopy. 

5.2.1 Results for the main stellar parameters 

Models with stellar masses within the 1 a log g range are 
typeset in boldface in Table [S] We found 6 models out of 
the 13 that fulfill this criterion. The best fit over all with 
Cr.jn.s. ~ 0.93 s has Mt = 0.65 Mq which corresponds with 
the expected mass from the spectroscopic log g value. Fig. [7] 
shows the 13 selected models in the Tcfi - log g plane. Closed 
and open circles denote the models within and out of the 1 a 
log g range, respectively. The spectroscopic solution with the 
corresponding uncertainties determined by Fontaine et al. 
(2001, 2003) is also given in Fig. [7] 

Assuming that the largest amplitude modes at 829, 866 
and 993 s are I — 1, we favour the models which give this 
solution. Within the 1 a logg range we find the 0.645, 0.650 
and 0.680 Mq models in accordance with this criterion. Tak- 
ing into account that TeH = 11 290 ± 200 K is expected also 
by spectroscopy, the 11 800 K solutions are over the 1 a limit 
and seem to be too hot. This means that the 0.645 Mq 
model would be our best choice. 

Table |6] shows the hydrogen layer masses of our models 
are between 10~* — 4 x 10~* M« . Restricting ourselves to the 
6 models selected by log g we find that 5 of them have Mh = 
10-''-6xlO-'^M.. We can give additional constraints on the 
hydrogen layer mass if we take into account only the three 
'favoured' models between 0.645 and 0.680 Mq. In this case 



Mh = 2.5x10" 



6.3xlO^^M.. 



5.2.2 Mode identification 

According to our selection criterion the models have mostly 
I = 1 solutions for the pulsation modes. Models with solely 



I = 1 modes exist only out of the 1 a log g range (M, = 0.725 
and 0.74 Mq). 

Considering the full range of best-fit models, we cannot 
uniquely assign an I value to any mode. At least one model 
has an Z = 2 mode solution for each period. We cannot 
provide unique assignments for the overtone numbers either, 
however, the range of stellar and hydrogen layer masses of 
our models constrains the k values. 

In summary we can say that the observed pattern of 
modes can be decribed with high accuracy assuming mostly 
I — 1 modes. Considering the 3 'favoured' models the modes 
at 619 and 778 s are always different I values. The 1250 s 
mode is always I — 2 for both these models. 



5.2.3 Varying the He layer mass 

In our second scan, we varied the mass of the He layer in 
the grid from 10~^ to 10~^ M* in 3 steps. This way we 
could check the sensitivity of our solutions to the thick- 
ness of the He layer. In this case we ran a more sophisti- 
cated version of the WDEC with core compo sition profiles 
based on stellar evolution calculations done bv lSalaris et al.l 
(| 19971 ). This grid covers the temperature and hydrogen layer 
mass range we used previously (Tcff = 10 800 — 11800K 
and Mh = 10"* - 10"* M,). We varied the stellar mass be- 
tween 0.60 - 0.69 Mq. Step sizes were as well 200 K (T^s), 
0.005 A/q (M.) and 10"°'^ M. (Mh). 

According to our first criterion we searched for models 
with at least 3 1 = 1 solutions and low tJr.m.a. values. We 
found two acceptable models at 0.61 and 0.635 Mq. Both 
models have the same helium layer masses (10"=*M.). How- 
ever, none of them have all of the large amplitude modes as 
1 = 1. The mode at 993 s is always 1 = 2, which would imply 
a much larger physical amplitude for this mode since it has 
the largest light amplitude. These two models have hydrogen 
layers with masses of 2.5x10"® and 10"^^,. These values 
are at the thin side of the range given by the previously 
selected models. Although we cannot rule out helium layer 
masses thinner than 10"^ M« , our favouring the three largest 
amplitude modes being I = 1 implies that Mho = 10"^ M« 
is the preferred value. 



5.2.4 Seismological parallax 

According to the luminosity of a selected model we can give 
an estimation on the distance and parallax of the star. Un- 
fortunately, there is no trigonometric parallax measurement 
for KUV 024644-3239. However, we can still predict what 
the parallax should be. 

The first step is to derive the bolometric magnitude 
{Mhoi) of the star: 



M, 



M 



- 2.5 log (L/Lq), (2) 

Qboi = + 4 75 (lAUenl Il973f) and = 
3.854 x 10 ergs/s (jSackmann. Boothrovd fc Kraenieil 

119931 ). Mboi can be related to the absolute visual magni- 
tude using the bolometric correction: My = Mboi — BC. 
Bergeron, Wesemael & Beauchamp (1995a) performed 
colour-index and magnitude calculations using DA and 
DB mo del grids. According to Table 1 in [Bergeron et al.l 
(|l995al ) BC = -0.441 and -0.611 at temperatures 11000 
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Figure 8. Period spacing diagrams of tlic 6 selected models by 
logg. Filled and open circles denote the model periods with 1 = 1 
or / = 2 solutions, respectively. Vertical dashed lines are drawn 
corresponding to the observed period values. Model periods found 
to be the closest of the observed ones are denoted with open 
squares. 

and 12 000 K, respectively. In view of the apparent visual 
magnitude jjIv = 16.07 and using the distance modulus 
formula we can derive the distance and the parallax of 
KUV 02464+3239. 

The 6 models selected by \ogg have logL/LQ = 
—2.56 to —2.71 where the most and least luminous mod- 
els are at 0.615 and 0.625 Mq, respectively. In these cases, 
the calculated seismological parallax values are between 
13.5 — 15.4 mas (milliarc second) with a 14.6 mas aver- 
age. Restricting ourselves to the 3 'favoured' models with 
log L/Lq = —2.6 to —2.66, the parallax angles are 14.8, 
14.2 and 14.7 mas for the 0.645, 0.650 and 0.680 Mq mod- 
els, respectively. These results imply a distance of ~70pc 
for KUV 02464-f 3239. 

5.3 Mode trapping 

KUV 02464-f 3239 was mentioned by iFontaine et"al] (120011 ) 
to be a photometric twin of another cool DAV, GD 154. 
In spite of the similarities between the shape of their light 
curves we find significant differences if we compare their 
pulsational properties. The Fourier spectrum of GD 154 
is dominated by three normal modes and their harmonics 



(|Pfeiffer et al.ll 19961). Subharmonics of the 1186 s mode were 
also detected by I Robinson et ahl (|l978h . In the case of KUV 
02464+3239 even with single-site observations we could de- 
termine 6 modes and using the whole light curve, only the 
first harmonic of the dominant mode was found. Clear pres- 
ence of further harmonics and subh armonics can be se en 
only by some parts of the light curve (|Bognar et al.|[2007^ . 

The small number of GD 154's eigenmodes compared 
to other cool DAV's implies a very efficient mode selection 
mechanism - possibly trapping of t he modes by a very thin 
(Mh ~ 10"^° M,) hydrogen layer (|PfeiflFer et al.lll996h . We 
can use our models on KUV 02464+3239 to test whether 
some of the modes are trapped. 

One possible way to find trapped modes is to construct 
period spacing diagrams and search for minima. Mode trap- 
ping causes departures from the uniform period spacings of 
consecutive k modes. Trapped m o des h ave lower kinetic en- 
ergies and as iBradlev fc Wingetl (|l99ll ) showed, minima in 
period spacing diagrams correspond to minima in kinetic 
energies. Accordingly we made a simple statistical fit on the 
6 models selected by log jr. We collected how often we find a 
period near or at a minimum of a model's corresponding pe- 
riod spacing diagram. We found that the largest amplitude 
modes (at 829, 866 and 993 s) frequently occur near or at 
period spacing minima and in these cases the I = 1 solution 
is preferred. The 1250 s mode has not been found at a min- 
imum in any of the cases. Fig. |8] shows the period spacing 
diagrams of the 6 selected models. Filled and open circles 
denote the model periods with Z = 1 or Z = 2 solutions, re- 
spectively. Vertical dashed lines are drawn corresponding to 
the observed period values. 



6 SUMMARY AND CONCLUSIONS 

We have presented the analyses of our observations on the 
DA variable KUV 02464+3239. 

Using analyses of data subsets we accepted 6 frequencies 
as normal modes of pulsation. With the analysis of the whole 
dataset we determined 7 additional frequencies as possible 
modes. Since remarkable amplitude variations (up to 30 - 
60 %) occured on a time scale of weeks, we performed tests 
for the possible sources of the amplitude variability. Neither 
beating of the 6 modes nor the effect of noise resulted in large 
enough variations. The test cases of closely spaced modes or 
rotational triplets revealed large amplitude variability, but 
only in unrealistic amplitude and phase relations. A real 
amplitude variation seem to be plausible. 

The best asteroseismological models have M* = 
0.645, 0.650 and 0.680 A/q and the mass of their hydrogen 
layer is between 2.5 x 10^^ — 6.3 x lO^'' M» . Although we can- 
not exclude thinner helium layers, A/ho = lO"'^ Mt models 
reproduce the observed properties of the star better. Aster- 
oseismological parallax values calculated by the luminosity 
of the 3 models are between 14.2 - 14.8 mas. Using period 
spacing diagrams we conclude that mode trapping can ex- 
plain the high amplitude modes, but is not required. 

The best way to discriminate between models would 
be to obtain observational constraints on the / values of at 
least the large amplitude modes. It may be possible by time 
dependent UV spectra or chromatic amplitudes but due to 
the faintness of the star and its complicated pulsation spec- 
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tra this would be a challenging task. Our results show the 
effectiveness of long, single-site observations. Multisite cam- 
paign(s) on KUV 02464+3239 would enable us to find addi- 
tional pulsation modes and possibly signs of rotational split- 
tings, and this would give further constraints on the stellar 
structure. 

Resolving the closely spaced modes expected in the case 
of KUV 02464+3239 we could exclude the beating of modes 
as cause of the observed amplitude variability. Our tests 
on amplitudes indicate this possibility. This would mean 
that we see real energy content variations of the individual 
modes. In this case the observed variability is presumably 
due to non- linear processes like mode coupling and/or inter- 
action between pulsation and convection. Finding the real 
physical explanation of amplitude variations is a key issue 
not only in the case of KUV 02464+3239, but also for other 
long-period white dwarf pulsators. 
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